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SUMMARY 

The methods of the calculus of variations are used to find functions of 
time for the thrust magnitude and direction that minimize the propellant con- 
sumption for a given Interplanetary transfer. The flight is assumed to be 
planar and to take place in an inverse -square gravitational field. Propul- 
sion is provided by a fixed thrust and specific impulse engine that can be 
turned on and off at will. With such an engine, the determination of the 
thrust -magnitude function reduces to finding the best points at which to begin 
and end a coast phase. 

Numerical results are presented for one-way trips from Earth to Mars and 
from Mars to Earth (both assumed in coplanar, circular orbits) for three 

values of initial acceleration (0.5x10“^, 1X10“^, and 5X10“^ Earth standard 
g r s) and five values of specific impulse (2000, 4000, 6000, 8000, and 10,000 
sec). For these engine parameters, the trip time Is varied between the 
continuous -propulsion time and the optimum time, which is analogous to the 
two-impulse, Hohmann transfer for impulsive thrust. These results are com- 
pared with the constant- jet -power, variable-thrust engine for the case of 
transfers to the best point on the Martian ellipse. 

Although the variational solution for round trips Is not considered, the 
one-way trip data of the type presented may be used to calculate equivalent 
results. A sample calculation is given to illustrate the procedure. 


INTRODUCTION 

The long propulsion periods associated with low-acceleration flight give 
rise to the problem of determining appropriate functions of time for engine 
thrust magnitude and direction in order to minimize the propellant consumed 
in accomplishing a given mission. As shown in references 1 to 4, such func- 
tions can be found by use of the indirect methods of the calculus of varia- 
tions. In formulating such a problem it is necessary to introduce constraints 


on the engine thrust and specific impulse consistent with the assumed mode of 
engine operation. Otherwise, as shown in reference 1, the variational solu- 
tion will yield infinite thrust and specific impulse. The purpose of this 
report is to present a variational solution with constraints specifically 
designed for early ion engines. This solution is then used to determine 
optimum trajectories for both the outbound and inbound heliocentric portions 
of low- thrust Mars missions. The resulting variational solution can also be 
applied to many other trajectory problems or to flight with chemical- and 
nuclear-rocket engines operating under similar contraints. 

In references 2 and 5 the thrust and specific impulse for ion engines were 
limited by constraining engine operation to constant beam power; the result- 
ing solution required large variations in thrust and specific impulse along 
the path. Early ion engines probably will be capable of operating efficiently 
over only limited ranges of thrust and specific impulse. Thus, although the 
variable-thrust trajectories of references 2 and 5 represent an upper bound 
on the capabilities of future ion engines, they are not applicable to early 
ion engines. The contraints assumed in the present analysis are that engine 
thrust and specific impulse (and, hence, beam power) are constant. Since it 
is assumed that the engines are capable of shutdown and restart, the tra- 
jectories may have coast phases. This problem is treated in reference 3 for 
the case of a constant planar gravitational field. Reference 5 presents a 
more comprehensive solution that includes a spherical inverse-square gravita- 
tional field and covers both the constant-thrust and variable -thrust cases. 

The solution of the present report is specifically for the constant -thrust - 
case in a two-dimensional inverse-square gravitational field and contains 
more detailed and extensive information about this case than has been pre- 
sented elsewhere. Emphasis is placed herein on a discussion of the nature of 
the optimizations and the effect of initial thrust -to-weight ratio and spe- 
cific impulse. Furthermore, typical results are given to illustrate the 
analysis of round-trip missions. 

The variational analysis results in a set of Euler-Lagrange differential 
equations that yield the optimum thrust direction and indicate the best times 
for initiating and terminating the coast phase. Numerical results were ob- 
tained by simultaneous numerical integration of the Euler-Lagrange equations 
with the equations of motion on an IBM 704 digital computer. As shown in 
reference 4, all the differential equations can be analytically integrated 
during the coast phase. The analytical solutions to the Euler-Lagrange 
equations for the present problem are presented in appendix B. The details 
of all numerical techniques employed here are given in reference 6. 

Sufficient numerical data are presented to indicate the nature of constant- 
thrust trajectories and to allow computation of some specific one-way and round- 
trip missions. Charts are presented for both Earth-Mars and Mars -Earth trans- 
fers that use an optimum travel angle. The data are presented as a function of 
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travel time and specific impulse for thrust-to-weight ratios of 0.5xl0~ 4 , lxlO~ 4 , 

and 5X10 ~ 4 . The const ant-thrust transfers are briefly compared with the 
variable -thrust transfers presented in reference 2. 

For round-trip missions , nonoptimum travel angles will generally be re- 
quired to meet planetary configurations. Thus, data for a wide range of 
travel angles are presented for a few representative values of thrust-to- 
weight ratio and specific impulse. These data are then used to illustrate a 
round-trip calculation. 


ANALYSIS 

The problem to be solved is to determine the vehicle thrust program and 
hence the trajectory that minimizes the amount of propellant consumed in ac- 
complishing a given mission. The engine is assumed to operate at constant 
thrust and specific impulse during powered flight, but it can be turned off 
for coast flight. For simplicity the vehicle is assumed to move in a two- 
dimensional trajectory in a central inverse-square gravitational field. 

The vehicle trajectory must satisfy the following equations of motion: 


u ? V i p 

= - In- 4- orr + — sin t 
r 2 m 

(la) 

2 U £0 j ^jP cos ^ 
03 ~ " r m r 

(lb) 

r = u 

(lc) 

Cp = 0) 

(Id) 

m = -3 

(le) 


(All symbols are defined in appendix A.) In order to satisfy the engine con- 
straints, the exhaust velocity Vh is treated as a constant throughout the 

cl 

flight, and the propellant flow rate [3 is given by 


0 = Pdes = const 

3 = 0 


for powered flight 
for coast flight 


( 2 ) 
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Formulation of Var ia tional Problem 


The variational problem to be solved is to determine f and g as func- 
tions of time, subject to the constraints given by equations (l) and (2), so that 
the amount of propellant consumed in accomplishing a given mission is minimized. 
This problem is a special case of the problem of Bolza as formulated in refer- 
ence 7 and is known in the calculus of variations as the problem of Mayer. The 
Multiplier Rule and associated corollaries developed in reference 7 will be used 
to solve this problem. In the problem of Mayer a function g of the initial 
and final conditions is minimized. To minimize the propellant consumption (or 
equivalently, the negative of final mass), g is given by 

g = -mf (3) 

The constraints given by equations (l) and (2) are introduced through the func- 
tion F defined by 


6 

F - Z Vi 
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where 


• u ? V i p 

> = u + it- - (xtt - — — sin if - 0 

J- r 2 m 

(5a) 

VjP 

fp = ro> + 2 ujd - — 1 cos = 0 

c m 

(5b) 

f 3 = r - u = 0 

(5c) 

f 4 = Cp - CD = 0 

(5d) 

f^ = m + g = 0 

(5e) 

f 6 - 3(e - 3 des ) = o 

(5f) 


and the A i (where i = 1, 2, . . . , 6) are the Lagrangian multipliers, which 
are unknown functions of time to be determined as part of the variational so- 
lution. 

Equations (5a) to (5d) are the equations of motion (eqs. (la) to (le)). 
Equation (5e) is the differential equation for the vehicle mass. Equation (5f) 
states the constraints on the propellant flow rate g given by equation (2); 
namely, that the propellant flow rate must be zero or the design value. At this 
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point it is convenient to introduce the following notation for the dependent 
variables of the problem: 


Zp(t) — U 

Z 5 (t) = m 

Zg(t) = 01 

z 6 (t) = i 

z 3 (t) = r 

z 7 (t) = p 

Z 4 (t) = 9 


Euler-Lagrange 

Equations 


As developed in reference 7, a necessary condition for g to be a minimum 
is that the Euler-Lagrange equation with respect to each dependent variable be 
satisfied. The Euler-Lagrange equations are given by 


d / dF \ _ 3F 
dt \dZjJ 


i — 1, 2, 


( 6 ) 


Applying equation (6) to equation (4) results in the Euler-Lagrange equations 
for the present problem: 


Ap = 2 odA2 ~ A 3 


(7a) 
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(7d) 

V,0 

A 5 = -*- (Ap sin t + A 2 cos 
mr 

*) 

(7e) 

^ v j p 

(Ap cos i|r - A 2 sin i) m : 

= 0 

(7f) 

V, . 

— ! a- (At sin + A 2 cos ty) 

m x 

- A 5 

(7g) 

' 6 = " 20 - 0 des 



The differential equations (7a) to (7e) must be integrated along with the equa- 
tions of motion (eqs. (la) to (le)) to determine the optimum trajectory. Equa- 
tion (7g) serves only to define Ag and has no further significance in the 
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present problem. Equation (7f) gives the value of the thrust angle along the 
optimum trajectory during powered flight (The value of the thrust angle is not 
required during coast flight). From equation (7f): 

A 1 

tan t = for p = P des (8) 

or 

A 1 

sin i|r = 

± V A i + A 1 

Ag 

cos = - 

± V A 1 + A l 

The uncertainty in the signs of equations (9a) and (9b) corresponds to an un- 
certainty of it in the value of the thrust angle i|r (i.e., i(r or t|r + it) . The 
choice of sign for equations (9a) and (9b) is developed later in the section 
Necessary Condition of Weierstrass. Integrating equation (7d) gives 

= const (10) 

Another integral of the Euler-Lagrange equations can be obtained, since 
the function F does not contain the independent variable, time, explicitly. 

From reference 7, 


(9a) 

(9b) 
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1=1 


In the present case, where the function F does not contain time explicitly, 
this reduces to 
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i=l 


or from equations (4) to (5f) 


($ - ■ -£r sln li + ( : 


V iP , 

2 un) - — * — cos I Ag 


0 


(12) 


uA-^ — cdA^ t |3Ac^ — 0 


( 13 ) 
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Equation (13) can be used to determine the value of one of the Lagrangian 
multipliers. In the present case it is convenient to determine A & from equa- 
tion (13) and eliminate the need for integrating equation (7e). Thus, 


C '( 

A 5 = ^ 


_y_ 

2 


co^r 


ViP 


m 


sin t 


\ / V,p \ 

1 Ap - I 2uo) - ^ cos I A 2 + uAg + o9\4 


for p ± 0 


(14) 


To evaluate A & for P = 0, it should be noted from equation (7e) that in this 
case A 5 = 0. Thus, A 5 is a constant during coast. Specifically, A 5 is equal 
to its value at the termination of the previous propulsive phase of the trajec- 
tory. The continuity of A 5 at the point of transition from powered flight to 

coast is established by applying the Weierstrass-Erdmann Corner Condition 
(ref. 7). This condition states that the SF/dZp must have equal right- and 
left-hand limits at such corners. For the present problem 


dF 

3u 




Sf 

3m 


a 5 


3f 


rA 2 



0 


3f 

3r 


dF 

^P 


0 


dF 


*4 


Therefore, at a corner such as the transition from powered flight to coast (and 
vice versa), the Lagrangian multipliers Ap to A^ have equal right- and left- 

hand limits. This implies that the constant of integration C is also continu- 
ous at the corners . 


Necessary Condition of Weierstrass 

Two uncertainties still remain in the system of Euler-Lagrange equations. 
The first is the choice of signs in equations (9a) and (9b), and the second is 
the value of p (i.e., P = P des or 0) along the path. Both of these uncer- 
tainties are resolved by applying the Necessary Condition of Weierstrass 
(ref. 7) . Following the development in reference 3, this condition states that 
for a minimizing trajectory E > 0 where 
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(15) 


/ * • * \ 

E = F(Zi, Zi) - F 


7 

(21, zi) - y 

i=l 


(z!-Zi)^>0 




The Zi's are the actual minimizing functions, and the Z*'s denote functions 
differing from the Zj/s by finite, admissible amounts. In the present problem 

only t and 0 are subject to such strong variations. Equation (15) may then 
be written as 


,pi 

m 


i r v i 

5 - P* 

1 m 


0 sin ^ + Ag cos i|r) - Ag - 0 * -*■ (Aj_ sin i|r” + Ag cos i|r*) - A 5 > 0 


5 > 


In the notation of reference 3, let 


V, 


k = (A-i sin i|r + Ap cos \|r) - Ac: 
m x ^ 0 


V, 


k* = .-U (Aq sin Tjr* + Ag cos ty*) - A 5 


(16) 

(17a) 
(17 b) 


Then 

Pk - p*k* >0 (18) 

Inequality (18) is now examined for the admissible variations in and p. 

P = P* = P(j_ es ; f * - In this case inequality (18) reduces to 


k - k* > 0 (19) 

or 


A-^ sin \|r + Ag cos \[r > sin \|/* + Ag cos \|r* ( 20 ) 

Since, for any value of ty, i|r* can only equal Tfr + jr , 

Aq sin Tjf + Ag cos i|r > 0 ( 21 ) 

Substituting equations (9a) and (9b) in inequality (2l) results in 

± + A| > 0 (22) 

Therefore, the plus sign is used in equations ( 9 a) and ( 9 b). 
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k(|3 - 3*) > 0 


(23) 


When k is negative, 3 must be less than 3* and thus 3 is zero. Simi- 
larly, when k is positive, 3 must be greater than 3* and thus 3 equals 
(3^ es . In summary , 


3 = 3a.es for k ^ 0 

3 = 0 for k < 0 


(24) 


The sign of k then determines the powered and coast phases of the flight, and 
the transition from one phase to another occurs at k = 0. 


The variational equations developed thus far are independent of the func- 
tion g; this is characteristic of the problem of Mayer, in which g is only a 
function of the initial and final conditions. The specific nature of g in- 
fluences only the boundary conditions that must be satisfied in integrating the 
equations of motion and the Euler-Lagrange equations. 


Boundary Conditions 

If some of the Z^s are not specified at either end of the trajectory, 

the free boundary values should be selected to minimize g. The boundary con- 
ditions that must be satisfied for the minimizing trajectory are given by the 
transversality condition (ref. 7). This condition states that for a minimizing 
trajectory 



By use of equations (4), (5a) to (5f), and (12), equation (25) becomes 

■£> 

C dt + Ai du + rAg do) + Aj dr + A^ dtp + Ag dm I + dg = 0 (26) 

l 0 

Further evaluation of the boundary conditions to be satisfied requires applica- 
tion of equation (26) to a specific problem. In this analysis the problem of 
minimizing the propellant consumed in transferring between circular, coplanar 
orbits is examined. (The problem of transfer between elliptic orbits Is 
discussed in appendix C). For this case g is given by equation (3), and the 
following initial and final boundary conditions are known: 
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Initial 


Final 


HlQ 
r 0 

= ynAo 

u 0 = ° 
cp 0 = 0 
t 0 = 0 

The final time and angle cpp may be specified or left free for further 

optimization of the final mass. As a consequence of the aforementioned boundary 
conditions (eqs. (27)), 

dciQ » d Xq = ckOQ =* duQ = dq>Q = dtQ = dr-^ = do^ = dup = 0 (28) 

and equation (26) reduces to 

C dtp + dcpp (At^ f ” -0 drip = 0 (29) 

where the variations dtp, dcpp, and dmp are independent. If tp, cpp, and mp 

are not specified, their variations in equation (29) are nonzero and their re- 
spective coefficients must be zero. 

Since mp is the quantity subject to optimization, it is not specified 
and dnip ^ 0. Thus, equation (29) gives as one boundary condition to be 
satisfied: 

A^ f = 1 (30) 

If the position in the final circular orbit is not specified dcp^ ^ 0, and 
equation (29) gives 

A 4 = 0 (31) 

Thus, a circular- orbit transfer with A 4 = 0 yields the trajectory possessing 

the optimum travel angle. If the final time is not specified dtp ^ 0, and 
equation (29) also gives 

C = 0 (32) 

Thus, a circular-orbit transfer with C = 0 yields the trajectory having the 
optimum travel time. 
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Calculation Procedure and Two-Point 


Boundary Value Problem 

The equations required for the determination of the minimizing trajectory 
are summarized for convenience In appendix D. The equations of motion (eqs. (la) 
to (le)) and the Euler-Lagrange equations (eqs. (7a) to (7c)) must be Integrated 
simultaneously to determine the minimizing trajectory. The value of the thrust 
angle Is given by equations (9a) and (9b), and the powered and coast phases are 
indicated by equation (24). These equations were programed for an IBM 704 com- 
puter and integrated numerically with a Runge-Kutta technique by use of a step- 
size control to limit the truncation error. 

For the problem of transferring between circular orbits the known initial 
conditions were given by equations (27). In order to proceed with the trajectory 
integration, values for the A i (0) must also be chosen. The values for the 

A i (0) must be selected to satisfy the desired final conditions given by equa- 
tions (27) and those obtained from the transversal! ty condition (eq. (29)); this 
results in what is commonly known as a two-point boundary value problem. 
Specifically, 


Unknown initial 
conditions 


Desired final 
conditions 


Xx(0) 

a 2 (o) 

a 3 (o) 

A^ = const 
C 


r(t f ) 
u(t f ) = 0 
a)(t f ) = f 

cp(t f ) 


(33) 


where for convenience C is selected as an initial condition to be determined 
rather than A^. 

To simplify the set of conditions (33), first note that the Euler-Lagrange 
equations (eqs. (7a) to (7g ) ) are homogeneous in A i . This means that the so- 
lution is independent of the initial magnitude of one of the A^ f s, which can be 
considered to act as only a scale factor for the other A i , s. In the present 
case, A3(0) was selected as the scale factor, and its Initial magnitude w^s 
taken as unity. After a solution is obtained, the requirement A 5 ^p — 1 can 
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then be satisfied by merely adjusting the scale factor Ag(0) , Thus Xj(0) 
and |> = 1 can be removed from the iteration. 

A further simplification vas made in the present analysis by not specifying 
cp^; hence, A^ vas removed from the iteration. The solutions for optimum cp^ 

vere obtained vith A^ = 0, and the solution for nonoptimum cpp were obtained by 

covering a range in A4. Finally, the desired t^ is obtained by terminating 
the trajectory integration at t = tp. The set of boundary conditions (33) then 
becomes the following: 


Unknown initial 
conditions 


Desired final 
conditions 


M°) 

^2(0) 

C 


r(t f ) = r f 

u(t f ) = 0 ► 

= Wj / r 3 

v ' f J 


(34) 


Obtaining a solution, then, requires that the values of three initial conditions 
be determined such that the three specified final conditions are satisfied. This 
problem vas solved by use of a three-variable Newton -Raphs on iteration scheme 
wherein the various required partial derivatives were evaluated by finite- 
difference methods. The iterations were terminated when residual errors in the 
end conditions were such that their removal would cost no more than a specified 
amount of propellant. The iteration scheme, as well as the criterion for termi- 
nating the iteration, is discussed in reference 6. 


RESULTS AND DISCUSSION 

The variational solution developed in the ANALYSIS has been used to investi- 
gate the heliocentric phases of Earth-Mars and Mars-Earth low-thrust transfer 
trajectories. In making preliminary analyses of low-thrust interplanetary 
missions, it is convenient to assume a series of two-body trajectories rather 
than treating the precise multibody problem (ref. 8 ). For example, for the 
Earth-Mars mission, the Earth escape spiral can be computed from consideration 
of Earth ! s gravitational field only. When the vehicle has traveled a sufficient 
distance from the Earth, the coordinates can be transferred to the sun; and 
during heliocentric flight, only the solar gravitational field can be considered. 
Again, when the vehicle has approached sufficiently close to Mars, the coordi- 
nates can be transferred to Mars; and during the descent spiral, only the Martian 
gravitational field is considered- 

The variational solution developed In the ANALYSIS can be applied to all 
three phases of such a journey. Only the heliocentric portion, however, has 
been investigated herein. Data for the. Earth escape and Mars descent spirals 
can be obtained from the charts of reference 8, which are based on a constant- 
tangential- thrust steering program. For constant -thrust spirals, this steering 
program has been shown to give near optimum results (ref. 2 ). 
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Boundary Conditions 


The boundary conditions imposed on the heliocentric trajectories presented 
in this report are based on transfer from geocentric to heliocentric coordinates 
when the vehicle has achieved escape energy with respect to the Earth and trans- 
fer from heliocentric to Martian coordinates when the vehicle is at escape energy 
with respect to Mars. Thus, the heliocentric trajectories begin with the vehicle 
in Earth's orbit about the sun and terminate with the vehicle in Mars' orbit 
about the sun. For simplicity, Earth and Mars are assumed to move in circular, 
coplanar orbits about the sun. Assuming a circular orbit for Mars can lead to 
significant errors in some cases. An elliptic orbit for Mars is treated in 
appendix C. The values in the following table were used for the orbits of Earth 
and Mars: 


Planet 

Radius , 

Polar angular 


r. 

velocity, 


m 

radians/sec 

Earth 

1.4950X10 11 

1.9910X10" 7 

Mars 

2 . 2779X10 11 

1.0586X10" 7 


These values are based on a solar gravitational constant |j = 1.3245X10 2 ® 

m'-’/sec^. The boundary conditions are satisfied to the extent that the residual 
errors result in an error in final mass of less than 0.1 percent of idq. 


Earth-Mars Transfers 

The transfer trajectories investigated herein are limited to direct trans- 
fers making less than 1 revolution about the sun. The resulting flights, in 
general, have either continuous propulsion or two propulsive phases with an 
intermediate coast phase. To indicate the nature of these transfers, the Earth- 
Mars transfers will be examined in detail for one set of engine parameters. 
Similar results would be obtained for Mars-Earth transfers. The reference set 
of engine parameters selected is an initial thrust- to-weight ratio of 1X10" 4 and 
a specific impulse of 8000 seconds. Later, data for both Earth-Mars and Mars- 
Earth transfers for a range of engine parameters are presented. 

Effect of travel time . - The effect of travel time on Earth-Mars transfer 
trajectories for the reference engine parameters is illustrated graphically in 
figure l(a). The flights shown have an optimum travel angle (corresponding to 
A 4 = 0). Tie quickest transfer, 181.6 days, is obtained with continuous pro- 
pulsion. As travel time is increased, the length of the coast phase increases; 
and at approximately 290 days, the trip having the optimum travel time (corre- 
sponding to C = 0) is obtained. This journey is analogous to the Hohmann 
transfer for impulsive thrust and has the longest intermediate coast phase. 
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(a) Trajectory profiles. 

Figure 1. - Effect of travel time on Earth-Mars transfers. Initial thrust- to-weight 
-4 

ratio, 1x10 ; engine specific impulse, 8000 seconds; Lagrangian multiplier X^, 0- 


Specification of longer travel times merely results in the optimum transfer with 
an additional coast phase either in Earth r s orbit before proceeding with the 
transier or in Mars 1 orbit after completing the transfer. Xf much longer travel 
times (corresponding to transfers making more than 1 revolution about the sun) 
are specified, transfers with more than one sequence of propulsion, coast, pro- 
pulsion could presumably be obtained. Transfers with such long travel times 
were not investigated in the present analysis. 

The variation of the thrust angle along the trajectories presented in 
figure l(a) is shown in figure l(b). The thrust angle for the C = 0 transfer 
(290 days) is nearly tangential (again analogous to the Hohmann transfer for im- 
pulsive thrust) . As travel time is decreased, increasingly larger radial com- 
ponents of thrust are introduced. The radial components are outward from the 
local horizontal for the near Earth phase of the trajectory and inward for the 
Mars phase. The thrust angle for the continuous -propulsion transfer is 
continuous, and the apparent discontinuity in figure l(b) is only a result of 
presenting the thrust angle from -jt to it rather from 0 to 2it. In general, 
dty/dt is positive during the propulsion phases of the flight. 

The effect of travel time on the final- to initial-mass ratio m^/niQ is 

shown in figure l(c) for the transfers with an optimum travel angle. As de- 
veloped in the ANALYSIS, the highest mass ratio is obtained for the transfer 
with C = 0. Initially, travel time can be reduced from the optimum with only a 
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small decrease in mass ratio. However, as the travel time approaches that for 
the continuous-propulsion case, the mass ratio decreases rapidly. 

The variation in the optimum travel angle with travel time is shown in 
figure 1(d). The planetary configuration required for an optimum-travel-angle 
transfer at a given travel time can be met once each synodic period. 

Nonop timum travel angles, A 4 £ 0 . - The transfers presented in figure 1 all 
possess an optimum travel angle. At each travel time it is also possible to 
make transfers with travel angles both larger and smaller than the optimum, at 
the expense of increased propellant consumption. Thus, for one-way missions it 
will generally be desirable to use transfers with an optimum travel angle. For 
round-trip missions of specified duration, however, it will generally be desira- 
ble to use nonoptimum travel angles for both the outbound and inbound segments of 
the mission; this is discussed later. 

The range of travel angles that can be achieved at a particular travel time 
for an Earth-Mars transfer is illustrated graphically in figure 2. The flights 
shown are those very near the maximum and minimum and the optimum travel angle 
for a travel time of 240 days. The limiting cases of maximum and minimum travel 
angle require continuous propulsion, and all the intervening travel angles are 
obtained with varying lengths of coast phase. The longest travel angle is 
achieved by initially moving inside Earth' s orbit before moving out to Mars, 

The shortest travel angle is obtained by moving out beyond the orbit of Mars in 
the latter phase of the trajectory and entering from higher radii. This phe- 
nomenon cannot be noted clearly because of the limited scale of figure 2 . 



Figure 2. - Comparison of Earth-Mars transfers with optimum and 
nonoptimum travel angles. Travel time, 240 days; initial 

thrust-to-weight ratio, lxl0“ 4 ; engine specific impulse, 

8000 seconds. 
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Travel angle, Tf, radians 

Figure 3. - Effect of travel angle on final- to initial-mass ratio 
for Earth-Mars transfers. Travel time, ?A0 days; initial thrust- 
to-veight ratio, lxlO" 4 ; engine specific impulse, 8000 seconds. 


The effect of travel angle on m^/r &q is presented in figure 3 for the 

reference set of engine parameters and a travel time of 240 days. Figure 3 
demonstrates the optimizing nature of A^ = 0 trajectories. Travel angles 

larger than the optimum value are obtained with positive values of A^, and 
travel angles shorter than the optimum are obtained with negative values of A^ 

Performance map for reference engine parameters . - From data similar to 
those presented in figure 3 for a range in travel time, a map of possible Earth 
Mars trajectories can be obtained. Such a performance map is presented in 
figure 4 for the reference engine parameters. For convenience, travel angle is 
presented as a function of travel time with the final- to initial-mass ratio as 
a parameter. The continuous -propuls ion flights form the outer boundary of the 
map. Flights inside the boundary have a coast phase, while those outside the 
boundary are not achievable with the reference engine parameters. Also, in the 
upper right portion of the map the flights have an initial as well as an inter- 
mediate coast phase. Similarly, in the lower right portion of the map flights 
occur with a final as well as an intermediate coast phase. The double- coast- 
phase trajectories are not new trajectories, but are trajectories obtained at 
shorter travel times and angles extended with an initial or final coast phase. 
The boundaries of these multiple coast regions have not been evaluated because 
they are characterized by k = 0 either at the start or end of the transfer, 
which is a very sensitive two-point boundary value problem that cannot be 
solved accurately without special machine program modifications. 




Travel time, t f , days 


Figure 4. - Effect of travel time on travel angle for various values of final- to initial -mass 
ratio for Earth-Mars transfers. Initial thrust -to -weight ratio, IxlO -4 ; engine specific 
impulse, 8000 seconds. 

Curves indicating the C = 0 and A4 = 0 trajectories are superimposed on 
the map; once again the optimum nature of these trajectories (as given by the 
transversality condition in the ANALYSIS) is apparent- The A4 = 0 curve con- 
nects the vertical tangents to the mass-ratio curves; this produces the optimum 
travel angle for each travel time- Similarly, the C = 0 curve connecting the 
horizontal tangents gives the optimum travel time for each travel angle. The 
highest mass ratio, 0.931, is obtained at the intersection of the A^_ = 0 and 
and the C = 0 curves. This is the 290-day transfer described in detail In 
figure 1; as mentioned previously, longer travel times were not investigated in 
the present analysis. 

Engine parameters . - The data presented thus far have all been for the 
reference set of engine parameters. The effect of thrust-to-weight ratio on 
Earth-Mars transfers with optimum travel angles is shown In figure 5 for a 
specific impulse of 8000 seconds. For each thrust-to-weight ratio, travel time 
is varied from the minimum time (continuous-propulsion flights) to the time 
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Figure 5. - Effect of travel time and thrust -to- weight ratio on final- to 
initial-mass ratio for Earth-Mars transfers. Engine specific impulse, 
8000 seconds; Lagrangian multiplier 0. 


maximizing the final- to initial-mass ratio (flights with C = 0) . The predomi- 
nant effect of thrust-to-weight ratio is the widening of the hand of achievable 
travel times as thrust -to-weight ratio is increased. In particular,, achieving 
short travel times requires relatively high thrust-to-weight ratios. Also_, at 
the long travel times, the optimum transfers (C = 0) approach the Hohmann trans- 
fer as thrust-to-weight ratio is increased. 

The effect of specific impulse on Earth-Mars transfer trajectories is shown 
in figure 6(a) for a thrust-to-weight ratio of 0.5X10" 4 . Again, the trajectories 
have optimum travel angles, and the travel time is varied from the continuous - 
propulsion case to that for a maximum mass ratio. The specific impulse, of 
course, has a major effect on the required mass ratio, but a change in specific 
impulse has only a minor effect on the trajectory used. For example, a change 
in specific impulse has almost no effect on the travel angle required for a 
transfer of given duration, as indicated by the nearly vertical constant travel- 
in g! e curves in figure 6(a) . A secondary effect of specific impulse is to widen 
the range of travel times available as the specific impulse is reduced. This is 
due to the greater mass expenditure and hence increased average acceleration at 
the lower specific impulses and is thus similar to the effect of thrust-to- 
weight ratio shown in figure 5. 
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Trajectory Charts for One-Way Missions 


Figures 6 and 7 summarize the trajectory data that have been obtained for 
Earth-Mars transfers. All the data are for transfers with an optimum travel 
angle , and the travel time varies from the minimum time for the continuous - 
propulsion case to the time for the maximum final- to initial -mass ratio. In 
figure 6, the final- to initial -mass ratio is presented as a function of travel 

time and specific impulse for thrust-to-weight ratios of 0.5X10~ 4 , 1x10 " 4 , and 

5X10 In figure 7, the apportionment of travel time into initial propulsion 
time, coast time, and final propulsion time is presented for the same range of 
engine parameters. Figures 8 and 9 present similar data for Mars-Earth tra- 
jectories with optimum transfer angles. Comparison of figures 6 and 7 with 
figures 8 and 9 indicates the similarity between Earth-Mars and Mars-Earth tra- 
jectories . 


Comparison with Variable-Thrust Trajectories 

In order to compare the results of this analysis with the results for the 
variable- thrust engine of reference 2, two modifications must be made. The first 
concerns finding the optimum I for a given jet power and travel time. The 
optimum I can be obtained as part of the variational solution as shown in ref- 

erence 5 or by simply constructing the constant jet power-to-initial mass ratio 
P/Mq envelope of curves at constant F/W q and I. The latter method has been 

used here and is Illustrated in figure 10. The second modification concerns the 
boundary conditions. For the comparison, Mars r orbit was assumed to be an el- 
lipse, of known eccentricity and semilatus rectum, coplanar to Earth 1 s circular 
orbit. Transfers to the Martian ellipse were made with an optimum travel angle 
to the best point on the ellipse for the given travel time. The transversality 
conditions for these end conditions are developed and discussed in detail in 
appendix C. 

For comparison, a nominal value of 0.04808 kilowatt per kilogram for P/Mq 
( which corresponds to (f/Wq)I = l) was assumed, and F/Wq and I were varied 

to construct the envelope curve. The resulting propellant fractions are shown 
in figure 10, where they are also compared with the variable-thrust engine re- 
sults. There is very little difference between the two curves at high values 
of t^, but the inherent advantage of the variable-thrust engine becomes more 

apparent as the travel time is reduced. This difference, at short times, be- 
comes about 5 percent of the initial mass. The significance of these differ- 
ences cannot fully be assessed until more is known about these engines. The 
propellant requirements of both engine types, however, are valuable as a guide 
In design and evaluation. 
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Figure 10. - Comparison of fixed-design engine and constant jet power engine for Earth-Mars 
transfers. Jet power-to-inltial mass ratio, 0.04808 kilowatt per kilogram; (F/tf 0 )l = 1. 


Round-Trip Calculations 


As pointed out previously, this report does not contain sufficient data for 
a comprehensive study of round-trip, Earth-Mars missions, A limited amount of 
nonoptimum- travel-angle data has been computed, however, and will be used to 
help illustrate a method for computing round trips that uses performance charts 
of the type shown in figure 4. In the calculation of round trips of given du- 
ration and waiting time, the travel time and travel angle for the outbound 
transfer may be selected at will. Once this has been done, however, the same 
variables for the inbound transfer must satisfy the rendezvous equations 


~ ^w ^out + ^in 

► 

^nPE “ 13 ^out + ^in - 


(35) 
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where 


t round-trip mission time 

t w waiting time 

t QU t outbound travel time 

tj_ n inbound travel time 

a>g Earth T s angular velocity 

Mars T angular velocity 
cp ou t outbound travel angle 

cp^ n inbound travel angle 

In view of the rendezvous equations it can be seen that optimum- travel -angle 
transfers can be used only for special combinations of t w and To opti- 

mize the round trip, for specified t m and t^ the pair of outbound parameters 
^out^ ^out^ mus ^ found that maximizes the mass at return to Earth. Such a 

problem is, in principle, within the scope of the indirect method of the calculus 
of variations. As such, however, it is beyond the scope of this report. 


To solve this same problem by use of performance maps such as the one pre- 
sented in figure 4, a selection is first made of F/W q, I, and 

From an appropriate performance map, the corresponding value of (m-p/mQ) can 

out 

be read. Thus the thrust-to-weight ratio at Mars arrival, (F/Wq) (iriQ/m^*) , will be 


known and can be used with data (such as those found In ref. 8) to find the mass 
ratio required to spiral into and out of an orbit about Mars. Thus, the value 
of (f/Wq) ^ can be computed and can be used together with t^ n and <p^ n as 

determined from equations (35) to obtain the necessary value of ^ from 

an appropriate inbound performance map. 


In order to minimize the number of performance charts required, the first 
choice of (t,cp) out can be made such that (f/WqK ^ has a value corresponding to 

that for one of the Inbound performance maps that is available. Furthermore, it 
should be noted from figure 4 that a variety of values for (t,cp) ^ can be se- 
lected at the same value of (m r /mrs) . Since all such cases will use the same 

1 u out 

inbound performance chart, it is possible to find the optimum value of (t,cp) ^ 


for a given value of (m^/imO by use of only two performance maps. This 

1 u out 

latter procedure is illustrated in figure 11, where all the values of (t,cp) . 


in 
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Final- to initial 
mass ratio, 
m f /m 0 

1 'th'.S t- 


I L I 


Continuous propulsion 


zxrr 


950-Day round trip, 

500-day waiting time 

( p /y 0 ut 

"210 220" 230 240f 

Travel time, tp, days 


Continuous propulsion 



Figure 11. - Effect of travel time on travel angle For various values of final- to initial-mass ratio for Mar3- 
Earth transfers. Initial thrust-to-weight ratio, 1.25xl0’ 4 ; engine specific impulse, 8000 seconds. 


along the dashed, curve were computed from equations (35) by use of t m = 950, 
t w = 500, and a series of outbound trips taken from figure 4 at 

(m^/m^)^^ = 0.875. The value of (m^/m^)^^ = 0.875 has been selected so that 

with (F/W q) = 1X10~ 4 , (f/Wq) = 1.25X10" 4 , and I = 8000 seconds (figs. 4 

and ll) , sufficient propellant was available to spiral into and out of an orbit 
about Mars at 1.1 X r^. From figure 11 it can be seen that the best trip of 

this type occurs when the dashed curve becomes tangent to one of the mass-ratio 
curves for the inbound trip. The details of the best case from figure 11 are 
given in the following table: 





Time, 

days 

Polar 

angle, 

radians 

Final- to 
initial-mass 
ratio, 

m f/ m 0 

Outbound 

209 

2.550 

0.875 

Spiral In 

36.9 

.337 

.954 

Park 

430.5 

3.938 

1 

Spiral out 

32.6 

.298 

.958 

Inbound 

241 

2.935 

.919 

Overall 

950 

! 10.058 

.735 


The case shown in the aforementioned table is the best possible for the chosen 

values of F/Wq, I, t m , t w , and (m^/niQ) . In order to find the optimum value 

out 

of (m-p/ir^) , the basic procedure illustrated by figure lx 1 have to be re- 

out 

peated for other values of (m f /m q) , and each new value would require a differ- 
ent inbound performance map. 014 ' 

One advantage of the method presented here is that it allows the computation 
of optimum round trips for many different combinations of t w and t m without 

any change in the number of performance maps required. Should it be necessary to 
investigate many different values of (F/Wq) and I, however, the number of 

charts required could become prohibitive. 


CONCLUDING REMARKS 

The calculus of variations has been used to derive differential equations 
that define planar, two-body, optimum trajectories for given boundary conditions. 
In the derivation, it is assumed that the gravitational field is of the inverse- 
square type and that the engine operates with fixed thrust and specific impulse 
but may be shut down and restarted at will. The resulting equations have been 
programed for an IBM 704 computer along with a three-variable Newton-Raphson 
iteration scheme to overcome the inherent two-point boundary value problem. This 
machine program has been successfully applied in an investigation of both Earth- 
Mars and Mars-Earth transfer trajectories with a range of values for Initial 
thrust-to-weight ratio and engine specific impulse. 

In order to compare results with those for continuously variable thrust as 
presented In reference 2, a special series of calculations has been made for 
transfers to an elliptical Martian orbit from a circular Earth orbit. The final 
comparison of these results shows that the two solutions differ very little at 
long travel times and that the variable -thrust solution has an advantage at short 
travel times that amounts to about 5 percent of the Initial mass for the case 
chosen. 
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Most of the data presented have the optimum travel angle at each travel time 
and are, therefore, most applicable to the study of one-way transfers. Suffi- 
cient nonoptimum-travel-angle data have been presented to illustrate a possible 
method for computing round-trip trajectories. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, September 4, 1962 
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APPENDIX A 


SYMBOLS 


C 

E 

e 

F 

F/W 0 

f 


g 

I 

J 


K 


X, 


£,uf 


k 


L 


m 

P/Mq 

P 

r 

t 


°m 


w 


u 

V, 


first integral to Euler -Lagrange equations 
Weierstrass excess function 
eccentricity 



initial thrust-to-weight ratio 
constraint equation 

function of initial and final conditions to be minimized 
engine specific impulse, sec 

functional to be minimized by variational methods 
constant of integration 

transversality functions for elliptic orbits 
function defining coast phases 

function associated with analytical solution for coast phases 
mass, kg 

jet power-to-initial mass ratio, kw/kg 

semilatus rectum, m 

radius , m 

time, days 

travel time, days 

mission time, days 

wait time, days 

radial velocity, m/sec 

engine exhaust velocity, 9-80665 I, m/sec 
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z 

a 

P 

Q 

A 

P 

9 

cp. f 

* 

CD 


dummy variable 

argument of pericenter, radians 
propellant flow rate, kg/sec 
true anomaly, radians 
Lagrangian multiplier 

3 / 2 

gravitational constant, m /sec 
polar angle, radians 
travel angle, radians 

thrust angle measured from normal to radius, radians 
polar angular velocity, radians/sec 


Subscripts : 
des design 


E 

f 

M 

m 

P 

w 


Earth 

final (t = t f ) 

Mars 

mission 

propellant 

waiting 

t = 0 


Superscripts : 

. d/dt 

' d/dcp 
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APPENDIX B 


SOLUTION OF EULEF -LAGRANGE EQUATIONS DURING COAST 

The purpose of this appendix is to derive analytical solutions for the 
following system of equations: 


Aj - 2cpA 2 A3 — 0 

(Bl) 

• . a 4 

- - ^2 + 2c P A l + — = 0 

(B2) 

t- (rAj + rCpA 2 + <pA 4 + C) 

(B3) 

r = r(cp) 2 - ^ 

v c 

(B4) 

rip = -2r(<p) 

(B5) 


which apply during the coast phase. Equations (B4) and (B5) have the well-known 
solutions : 


r 


P 

1 + e cos (9 




(B6) 


and 



(B7) 


Changing the independent variable in equations (Bl) and (B2) from t to cp 
leads to 



- 2cpX£ + = 0 


(B8) 


and 



Ag + 2<pA 1 + 



(B9) 


where A^ = dA^/dcp., and so forth. Substitution of equation (B3) into ( B8) 
results in: 

AnCpr - (2cpr + rcp)A 2 - rA 4 - C - q>A 4 = 0 


32 



From equation (B5), 


2q)r + rqj = 0 


Therefore, 


T 




(BIO) 


By use of relations (B6) and (B7), 


equation (BIO) can be integrated to give 



L sin 0 
qe 



cos 0 + K-^ 


sin 0 


where 


L = 



de 

sin2©(l + e cos 0)2 


(Bll) 



sin 0 

(l + e cos 0) 2 


(l + 2e 2 )sin 0 
(1 - e 2 )(l + e cos 0) 


cot Q 

(l + e cos 6 )^ 


for 0 < e < 1 



G = cp - a and Kq is a constant of integration. In evaluating L , only the 
elliptical case has been presented because of its frequent occurrence. Substi- 
tuting equation (Bll) into equation (B9) and integrating lead to 



Q)^ + cot G 


j + ^4 ( f + cos e) 


fZ K-^l + e cos 0) 2 'i 

x Sln 0 Vq e K 2 f (B12) 

where Kg is a constant of integration. Equations (B6), (B7), (Bll), (B12), 
and (B3) can then be used to determine A3. 


Once the initial values of A^> Ag, 0, p, and e are known, the constants 
Kq and Kg can be evaluated with equations (Bll) and (B12) . The problem that 
still remains is to solve for the A^_ at the second k = 0 point, the termi- 
nation of coast. Since the A^ are given as functions of 0, the value of 0 
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at k = 0 must be determined. Because of the transcendental nature of the 
equations, some iterative solution must be used to calculate 0 at the termina- 
tion of coast. A method based on finite differences, which converges rapidly to 
0 at k = 0, has been successfully applied. The subsequent power phase is 
numerically integrated by the methods previously discussed. 
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APPENDIX C 


ELLIPTIC ORBIT BOUNDARY CONDITIONS 


Transfers to an elliptic orbit are in principle no different from transfers 
to a circular orbit; however, certain parameters may be more convenient to use 
than others. In this analysis, the eccentricity, semilatus rectum, and perigee 
angle are used to describe the shape, size, and orientation of the ellipse with 
respect to the initial circular orbit. The true anamoly is used to identify the 
position of the vehicle in the final orbit. It should be noted here that the 
orientation of the ellipse gives rise to another degree of freedom - a circle has 
no orientation. The transformation of variables can be made through the use of 
the following equations: 


u =* e sin 6 


03 = (l + e cos O)^ 


> 


r = 


(Cl) 


1 + e cos 6 
cp = 6 + a 

Substitution of the total derivatives of equations (Cl) into the transversality 
expression (eq. (26)) gives 


( 


A 5 ^ 


dJ = -dmp + A5 dm + ( — Ai + 


2r^co cos 8 , cos 8 

— Ao 


3 ) 


A 3 )de 




^1 2r^eco sin 8 


tan 8 


A 2 + 


r^e sin 0A 3 


+ A^JdS + A^ da + C dt 


(C2) 


or, for the case of fixed initial conditions, 

dJ = |j A 5 - l) dm + ){ de + £ dp + d0 + A4 da + C dt^J = 0 


(C3) 


where iK, £, and *Jl are coefficients of de, dp, and d0 as defined by equa- 
tion (C2) . As previously discussed, a transfer is optimum with respect to a 


35 


particular variable if its coefficient is zero at the boundary. Consider, for 
example, transfers used in the comparison with the variable -thrust engine of 
reference 2. Here, each transfer was made with e f , pp, and tp fixed; as pre- 
viously discussed, A 5 p = 1 was satisfied by scaling. In addition, Jl = 0 and 
= 0 were also satisfied. These transfers then have the optimum orientation a 

of the Martian ellipse (e,p), and the vehicle arrives at the optimum point on the 
final ellipse 0 f in the given time t f . Because of the relation among orien- 
tation, true anomaly, and travel angle, the optimum orientation also implies op- 
timum travel angle for some given The two-point boundary value problem for 

this case is, therefore, 


At t = t r 


*1 

Ao 


At t * tf ^ 

e = = 0.093369 

P - Vy[ ** 2.2581X10“^ 

J(- 0 


(C4) 


Another type of transfer that can be made is the same as that just described 
except that the point of arrival is nonoptimum 6. The optimum point of arrival 
jl~ 0 is included in this class of transfers. This class of transfers also 
points out that there are two points on the final ellipse that have ^ 0 - 
maximum and minimum mp. Therefore, an additional test should be made to 

distinguish the two cases. Also, some points on the Martian ellipse are not 
achievable with all sets of engine parameters and are defined by the continuous - 
propulsion limits of the engine. In this case, the minimum mp occurs at the 
continuous-propulsion points but is not defined by Jl = 0. These characteristics 
of the elliptic orbit boundary are shown in figure' 12, where both Jt and mp/mQ 

are given as functions of 9p for an Earth-Mars transfer with F/Wq =» 1X10" 4 , 

I = 10,000 seconds and t f = 200 days. 
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Final- to initial-mass ratio, m^/ra 0 



Trans versality function for 
elliptic orbits, Jt 





APPENDIX D 


EQUATIONS FOR MINIMIZING TRAJECTORY 

The following equations are required for the determination of the minimizing 
trajectory: 


Equations of motion: 


Euler -Lagrange equations: 


u = - -^r + co^r + — sin ^ 
r 2 m 


- 2ucd , cos ^ 

cd =* - + — - 

r m r 


Cp = CD 


m = -3 


A^ — 2cdA^ “ A^ 


* U ' '4: 

A g = -2 cd/v]_ + “ Ag - 


^ + do 2 J Aj_ + <dA 2 


where 


sin ^ and cos t are given by 


= const 


sin f 


+ J-h\ + a| 


cos ^ 


^1 + A l 
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0 is given by 


P = Pdes for k ^ 0 


0 = 0 


for k < 0 


V 4 

k = — ^ (Aj sin t + ^2 cos “ ^5 


~(^2 ~ aj2p ~ Sln ~ CQS + ^5 + °^4 


for 3^0 
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